Analytic treatment of a trading market model 
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We mathematically analyze a simple market model where trading at each point in time involves 
only two agents with the sum of their money being conserved and with neither parties resulting 
with negative money after the interaction process. The exchange involves random re-distribution 
among the two players of a fixed fraction of their total money. We obtain a simple integral nonlinear 
equation for the money distribution. We find that the zero savings and finite savings cases belong to 
different universality classes. While the zero savings case can be solved analytically, the finite savings 
solution is obtained by numerically solving the integral equation. We find remarkable agreement 
with results obtained by other researchers using sophisticated numerical techniques 0]. 



The distribution of wealth in a nation among its people 
has been of considerable interest over the ages. If one suc- 
ceeds in mathematically modeling the distribution it has 
serious implications. One of the empirical laws known is 
that due to Pareto 2] which states that the higher income 
group distribution decays like a power law with exponent 
between 2 and 3. On the other hand the lower income 
group distribution is Gibb's like. The Gibb's law has 
been numerically demonstrated to be obtainable when 
trading between two agents is completely random and 
docs not involve any savings |^ ^ |^ . Also the finite sav- 
ings case has been studied numerically by Chakraborti 
and Chakrabarti |^]. However an analytic treatment of 
these models is lacking and hence the present work is 
aimed at meeting this need. 

In this paper we derive analytically an integral equa- 
tion for the probability distribution function of the 
money in the system. The model describes trading 
between two agents i and j each with money mi{t) 
and rrijlt) at instant t. Each agent saves a fraction A 
of his or her money and the rest is traded. Trading 
involves a random re-distribution of the money avail- 
able for trading [(1 — X){mi{t) + mj{t))] among the two 
agents i and j. After trading , i.e., at instant t + 1, 
mi(t -|- 1) = •mi{t) -f Am and Tnj{t -f 1) = rrijlt) — Am 
where Am = [eijrii -\- rrij) — mi(t)](l — A) with < A < 1 
and e being a random number between and 1. 

We will now derive the equilibrium distribution func- 
tion f(y)dy which gives the probability of an agent having 
money between y and y + dy. We assume that, irrespec- 
tive of what the starting point is, the system evolves to 
the equilibrium distribution after sufficiently long time. 
We will now consider interactions after the system has 
attained steady state. The joint probability that, before 
interaction, money of i lies between x and x + dx and 
money of j lies between z and z + dz is f{x)dxf{z)dz. 
Since the total money is conserved in the interaction, we 
let L ^ X + z and analyze in terms of L. Then the joint 
probability becomes f {x)dx f {L — x)dL. We will now gen- 
erate equilibrium distribution after interaction by noting 
that at steady state the distribution is the same before 
and after interaction. Probability that L is distributed 



to give money of i between y and y -I- dy is 



dy 



(l-A)L 



f{x)dxf{L — x)dL, 



(1) 



with x\ < y < x\ + {1 — \)L. Thus we see that x < y/X 
and X > [y — (1 — A)iy]/A. Actually x should also sat- 
isfy the constraint < x < L because the agents can- 
not have negative money. Thus the upper limit on x is 
min{i,y/A} (i.e., minimum of L and y/X) and the lower 
limit is max{0, [y — (1 — A)i]/A}. Now, we know that the 
total money L has to be greater than y so that the agents 
have non-negative money. Thus we get the following dis- 
tribution function for the money of i to lie between y and 
y + dy 



f{y)dy = dy 



dL 



y (1 - A)L 

m'm[L,y / A] 



dxf{x)f{L-x). (2) 



max[0,{a-(l-A)L}/A] 

It is of interest to note that when A = 0, the upper and 
lower limits of the x integration become and L and we 
get 



/(y)dy ^dy J ^ ^ 



dxf{x)f{L - x). 



(3) 



We will now provide an analytic solution for the above 
zero savings case. We first note that the double derivative 
with respect to y of the above Eq. Q yields 



f'{y) + yf"{y)^^fiy)f{0) 



f{x)f{y-x)dx. (4) 



For small y we assume that the function /(y) and its first 
and second derivatives are well behaved. Then as y — > 
we get 



/'(y) « -/(y)/(0). 



(5) 



Then the solution for small y, after using the constraint 
ir fiy)dy = is given by 



/(y) « f{0)exp[~yf{0)]. 



(6) 
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Next we also note that the above function (in Eq. 
also solves the parent Eq. ||2J) exactly. It is of interest to 
note that there are similarities between this distribution 
and the classical Maxwell— Bolt zmann distribution with 
the maximum corresponding to y = 0. 

We will now make a few observations regarding f{y). 
On examination of Eq. (j^J, assuming /(y) — > as y — > oo 
and that /(y) is well behaved as y ^ 0, it is clear that 
limy^o f{v) 0. Physically this makes sense because if 
everyone has non-zero savings, then a person with zero 
money will, due to interactions, tend to non-zero money 
faster than returning to zero money. A person with zero 
money after a single interaction has probability 1 of get- 
ting non-zero money. Also, once a person has non-zero 
money it takes an infinite number of interactions for that 
person to end up with zero money. Thus the A = case 
and A > case belong to different universality classes. 
Moreover A = 1 also belongs to a different universality 
because it is a static situation. Also it can be shown that 
for 0<A<lasi/^C)o the function decays exponen- 
tially. 

We have solved the integral equation given by Eq. Q 
for the non-trivial case of f{y) ^ <5(2/)- Our approach is 
based on solving for f(y) iteratively by starting with a 
trial function. To achieve convergence, we use a novel 
approach, wherein as soon as /(y) is calculated at any 
y it is immediately used to evaluate /(y) at other values 
of y. Thus we do not need to wait till /(y) is evaluated 
completely for all values of y to use it to get the next 
approximation for /(y) as one might naively expect. De- 
tails of this iterative procedure and a critical analysis will 
be presented elsewhere. We give here results (see Fig. 
1) for integrations done with rectangular boxes of width 
0.01 and upper limit of L = 10. All the curves are nor- 
malized and correspond to average money being unity. 
As A increases, for A > 0, the peak value of the bell- 



shaped curves increases and shifts to the right and also 
the width decreases. It is of interest to note the rather 
striking similarity between our plots and those reported 
by A. Chatterjee, B. K. Chakrabarti, and S.S. Manna 
1] . It is very gratifying to note that by using completely 
different approaches one gets very similar results which 
lends confidence to these approaches. 

In future we hope to extend our treatment to obtain 
a model that is consistent with Pareto law. Numerical 
advances in obtaining the Pareto law have been reported 
inRef. 6]. 

The authors are very grateful to B. K. Chakrabarti for 
inspiring and useful discussions. One of the authors (S. 
Y.) also thanks R. Ramakumar for support and helpful 
discussions. 




FIG. 1: Plot of the money probability distribution function for 
various savings values (A = 0.0,0.1,0.5,0.9). The average money 
per person is set to unity. 



[1] A. Chatterjee, B. K. Chakrabarti, and S. S. Manna in 
these proceedings. 

[2] V. Pareto, Le Cours d'Economique Politique, Lausanne & 
Paris, (1987). 

[3] B. K. Chakrabarti and S. Marjit, Ind. J. Phys. B69, 681 
(1995). 



[4] A. Chakraborti and B. K. Chakrabarti, Eur. Phys. J. BIT, 
167 (2000). 

[5] A. A. Dragulescu and V. M. Yakovenko, Eur. Phys. J. 
BIT, 723 (2000). 

[6] A. Chatterjee, B. K. Chakrabarti, and S. S. Manna, 
|cond-mat/0301289| (2003) . 



